Multi-parameter generalization of nonextensive statistical mechanics 
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We show that the stochastic interpretation of Tsalhs' thermostatistics given recently by Beck 
[Phys. Rev. Lett 87, 180601 (2001)] leads naturally to a multi-parameter generalization. The 
resulting class of distributions is able to fit experimental results which cannot be reproduced within 
the Boltzmann's or Tsallis' formalism. 
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Nonextensive statistical mechanics (NESM) intro- 
duced by Tsallis has gained a considerable interest 
in several fields of physics because of its capability to de- 
scribe a wealth of disparate phenomena (from anomalous 
diffusion, to turbulent systems, to astrophysical systems, 
etc ...) within a single formalism, generalization of the 
standard statistical-mechanical one with the addition of 
the single free parameter (entropic index) q. Recently it 
has been shown how to relate q with the internal micro- 
scopic properties of the system under consideration. This 
has been done by Wilk and Wlodarczyk Q: they have 
shown that, when g > 1, the NESM canonical distribu- 
tion pq[H,(3f)) for the system with Hamiltonian H can 
be written as an average of the usual Boltzmann-Gibbs 
factor over the inverse temperature (3: 



(1) 



where /g(/3, /3o) is a weight function whose meaning is 
that of a probability distribution function for P which 
is, therefore, no longer a fixed parameter; instead, the 
macroscopically visible value is just its average value Pq. 
Fluctuations in f3 are related to coherent fluctuations ex- 
isting in small parts of the system in respect to the whole 
system, due to the existence of long range correlations. 
Recently, Beck ^ has been able to give an interpretation 
of the fluctuating /3 as a function of stochastically varying 
microscopic variables. In order to recover Tsallis' results. 
Beck was forced to impose some constraints over /? or, 
equivalently, the microscopic dynamics of the system. In 
this paper we show that, following Beck's approach but 
relaxing these constraints, we are able to derive an entire 
new class of distributions, which reduce to Tsallis' dis- 
tribution under suitable limits. We will show that some 
members of this class are able to reproduce experimen- 
tal results which would be outside the reach of Tsallis' 
formalism. 

To start with, we quote the same example used in 
Beck's paper: let us set H = and suppose that 

the generalized velocity u satisfies the Langevin equation 



with L(t) Gaussian white noise of unit amplitude, a 
strength of the noise, 7 friction coefficient. This is the 
Brownian particle problem [Q. For this case, it can be 
shown that the temperature is related to the micro- 
scopic parameters 7, a by 



(3) 



Beck shows that Tsallis's distribution can be recovered if 
f3 is characterized by a distribution with n degrees of 
freedom [Hi 
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where /g(/3, Po) — fn{P, Po) provided that q = 1 + 2/ {n + 
1). Such a distribution arises if P can be written as a 
sum of stochastic variables: 



(5) 
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with < Xi >= and < Xf >= Po/n, so that < P >= Po 
and < P^ > - < P"^ >= Pq (2/n). The distribution is 
a common distribution, occurring in many physical prob- 
lems, and is central in the problem of estimating param- 
eters from data Q. 

Some points are worth stressing at this stage: 

(i) The macroscopic parameter P is written in terms of 
other parameters more directly related to the microscop- 
ical dynamics of the system at hand, just as in Eq. (|^). 
We just mention another example: in the study of fully 
developed turbulence, where w is a local velocity differ- 
ence, P = (er)"^, with e spatially averaged energy dissi- 
pation rate and r typical time for the energy transfer. 

(ii) It is obvious that, if /3 is a stochastic variable, a 
fortiori the microscopic quantities 7, cr, . . . must also be 
stochastic variables, therefore characterized by their own 
probability distribution functions (PDFs). (iii) Relations 
of the kind (^ impose severe constraints upon the PDFs 
of the microscopic variables. For example, to recover Eq. 
(0) starting from Eq. (0) there is the trivial choice: 7 
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X^-distributed and cr^ a constant; it is difficult (and per- variables; after integration over we get 



haps impossible) to devise other distributions which lead 
to Eq. (1). 

The main idea of this paper is that if /3 is a function of 
some more fundamental stochastic control variables, then 
the by far more logical path is the following: to guess sta- 
tistical distributions for the microscopic quantities and, 
from them, to work out the corresponding distribution 
for p. Since /? may have infinite functional dependences 
from microscopic variables, we can expect the PDF of 
/3 to have a large range of analytical forms, depending 
on a large number of parameters [we expect as many of 
them as the number of microscopic variables that control 

/3 = /3(r,^,...)]. 

Some simple rules, however, still allow to drastically re- 
duce the class of likely distributions. First, although the 
PDF for each of the variables Y,Z,... , may be arbi- 
trary, the same reasoning of Eqns still holds: that is, 
the distribution for each variable is a very convenient 
choice. For example, the distribution can tranform 
into a delta distribution, thus allowing for well deter- 
ministic, non-stochastic quantities in the limit n oo. 
Hence, we will suppose all the stochastic variables to be 
X^-distributed, possibly with different degrees of free- 
dom. In second place, a simplicity principle suggests 
that the most frequently occurring cases should be those 
where P is some simple combination of a small number 
of variables. Some examples are given in the above ex- 
pressions (e.g., Eq. 1^). The simplest function of all is the 
sum of stochastic variables: f3 — Y + Z + ... . However, 
with the previous choice for the PDFs oi Y, Z, . . . , it is 
possible to show that it is a trivial case, since it reduces 
to a distribution The next nontrivial cases, thus, 
are those involving products and ratios of one or two con- 
trol variables: Y x Z, Y/Z , 1/{YZ), .... 
Our aim now is to compute a few examples of PDFs of /3 
and to compare the results with the Tsallis' formalism. 
We will do the computation for the case of f3 ratio of 
two stochastic variables: (3 = Y/Z. This is particularly 
convenient since: (i) it generalizes the example given by 
Beck (Eq. ||) ; (ii) it is a particular case of /3 = 1/ (er) , 
when Y and either e or r are constants. 
The probability distribution function for the two x^ 
dependent variables Y, Z of degree n, m respectively, is 
given by 
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(r is the factorial function r(z) = *dt). We set 

P = Y/Z, Po = Yq /Zq and regard P and Z as independent 
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which is known as F distribution in statistics. This is the 
main result of the work, since the statistical properties of 
the system are determined through the two-parameters 
canonical distribution, generalization of Eq. (|l|): 
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The main feature of Eq. (0) is that the exponential term 
of Eq. (^) has disappeared, replaced by a power-law 
term. One should expect this term to depress high-energy 
tails in Eq. (^) . In order to have an insight on the trends 
of Eq. (|^), let us consider some interesting limits. First 
of all, we observe that, in the limit m ~f oo, 
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We recover the distribution (Eq. U) since, in the 
limit of infinite degrees of freedom, the distribution for Z 
shrinks to a delta distribution, so we are actually dealing 
with just one stochastic variable, Y. It is completely new 
the limit n — > cxd (that is, we are computing the PDF of 
the variable l/Z), for which we get: 
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In order to give a visual insight, we plot in Fig. |l| some 
examples of these distributions. The qualitative shape 
is rather similar. The occupations factors are computed 
through eq. (^). We give explicit expressions for the 
cases corresponding to the two limits n — > cx), m — > oo: 



{H,Po) 
(H,po) 



{2mPoH)'^ 



2f-ir(^) 



(11a) 

^/2mPaH] (lib) 



where K is the modified Bessel function of order to/2 and 
erfc is the complementary error function: erfc = 1 - erf. 
The general case of arbitrary n, m can be explicitly writ- 
ten down, but it is not revealing since it involves complex 
combinations of hypergeometric function, difficult to vi- 
sualize. We plot in Fig. ^ the standard Boltzmann-Gibbs 
factor together with the curves (pi]). In general, the new 
distributions are characterized by tails intermediate be- 
tween Boltzmann's and Tsallis' statistics. We can obtain 
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FIG. 1: Probaility distribution fn,m{l3) from Eq. (|), with 
/3o = 1- Solid line, /a, 3; dotted line, /s.oo; dashed line, /oo,3. 




FIG. 2: Generalized canonical distributions pn.m as a 
function of the scaled energy PoH. Solid line, n — » 
00, m — > CX3 [this yields the usu al Bo ltzmann-Gibbs (BG) case 
exp(— /SoJy)]; dotted line, Eq. ( |lla[ ) with n = 1 , correspond- 
ingto the Tsallis distribution with q = 2); dashed line, Eq. 
(Lib) with m = 1 \ dotted-dashed line, n = 2,m = 1. 



the probability distribution Pn,m{u) for the generalized 
velocity u once an explicit form for H = H{u) is given. 
By assuming the usual form H = 



-Pn,oo (^) 



Po r(t) 
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Notice that the function K yields a typical exp(— cff 2 ) or 
exp(— c'|m|) dependence. Such a law cannot be recovered 
within the Tsallis' formalism, which predicts power-law 
dependences. Therefore we take it as a signature of this 
new class of functions. It may be of interest to notice 
that the dependence on \u\ comes from the variable at 
the denominator of while the numerator provides a 
dependence on u^. In the general case, both \u\ and 
terms do appear. 

The question arises obviously if such distributions do ex- 
ist in nature. We are interested in fluctuations of some 
quantity: for independent fluctuations, the Central Limit 
Theorem predicts a Gaussian PDF. If departures from 
Gaussianity are described in terms of Tsallis' statistics, 
only PDFs with power-law asymptotics may be included. 
On the basis of what told before, we must look for PDFs 
with exponential tails. Actually, in literature are pre- 
sented several examples of quantities whose PDFs are 
(at least on some ranges) exponential. We briefly men- 
tion the numerical computation of the velocity distribu- 
tion function solution of the Enskog-Boltzmann equation 
for a granular gas [Q; other hints come from calcula- 
tions of the large-scale probability density distribution 
in astrophysics and from the numerical simulation of 
stresses in sheared granular materials A field where 
several well documented examples can be found is the 
study of turbulence in fluids. We refer in particular to 
papers |jl^, ^ |l^, ^). The quantity we are interested 
in here is the PDF of the velocity difference between two 
spatial points. It is found both experimentally and nu- 
merically that this quantity shows an exponential tail. 
In particular, in paper pO| | the departure from a Gaus- 
sian form is interpreted within a formalism very close to 
ours, where the average (^ is done using their equivalent 
of fn,m{P,Po) given by a log-normal function [see their 
equations (3.1-3.4)]. The paper [Q, furthermore, shows 
that the tails of this PDF can smoothly vary between 
the Cauchy form (which is a particular kind of Tsallis' 
distribution) to Gaussian form passing through the expo- 
nential form, by varying a few control parameters. This 
is strikingly reminiscent of varying n, m parameters in 
our formalism. 

In more detail, we can quote two experimental stud- 
ies from fusion plasma physics: in the first paper |l^ 
it is presented a study of the density fluctuations ex- 
isting in a thermonuclear fusion device. The time be- 
havior of the electron density Ue close to the bound- 
ary of the device was measured with high sampling fre- 
quency, thus allowing to compute the PDF of the fluc- 
tuation fie = Ug— < He >. It was found that the curve 
is highly asymmetrical, with the negative wing approxi- 
mately gaussian, and the positive one nearly exponential. 
In Fig. H we fit the experimental data with both Tsallis' 
and our curve, showing that the former curve cannot fit 
the tail of the experimental distribution. Rather closely 
related, we mention a second paper, dealing with a statis- 
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FIG. 3: Probability distribution P{hs) of the electronic den- 
sity fluctuations he. Broken line, experimental data from ref. 
mM (only the side of positive fluctuations is shown); dotted 
line, best fit using Tsallis' distribution (^^; dashed line, best 
fit with curve ( |l3| ) and m = 1. 



was recovered by fitting them with stretched exponen- 
tials: P{X) w exp(— The parameter a is func- 
tion of the time scale, varying between 1 (exponential 
distribution) and 2 (Gaussian distribution). In Fig. ^ 
the case closest to an exponential is shown. 

We think we have given in this work constructive ev- 
idence of the existence of generalized nonextensive dis- 
tributions. The very simple PDFs we have computed, 
seemingly gave us the tools to describe complicated phe- 
nomena. 

A crucial point is the choice of the microscopic variables, 
since one could always choose varying definitions for them 
so as to identify several different cases within the same 
classes of functions. Therefore, work in this direction 
should: (i) either show that trivial redefinitions of vari- 
ables are not important for the final result, (ii) or find 
that some sets of variables are preferred with respect to 
all the others. 

We gratefully acknowledge G. Antar for providing us 
with the experimental data of Fig. H, E. Martines for the 
data used in Fig. ^ 
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FIG. 4: Probability distribution P{V) of the electrostatic 
potential fluctuations V. Circles, experimental data from ref. 
[p^ ; solid line, best fit with curve ([l^ ) and m — 1. 

tical analysis of electrostatic potential fluctuations, still 
in the edge of a plasma A wavelet analysis of the 

data allowed there to compute PDFs as function of the 
time scale of the fluctuations. A scaling law for PDFs 
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